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ABSTRACT 

We present a direct detection of the gravitational lens that caused the microlensing event MACHO- 
95-BLG-37. This is the first fully resolved microlensing system involving a source in the Galactic 
bulge, and the second such system in general. The lens and source are clearly resolved in images 
taken with the High Resolution Channel of the Advanced Camera for Surveys on board the Hubble 
Space Telescope (HST) ~ 9 years after the microlensing event. The presently available data are not 
sufficient for the final, unambiguous identification of the gravitational lens and the microlensed source. 
While the light curve models combined with the high resolution photometry for individual objects 
indicate that the source is red and the lens is blue, the color-magnitude diagram for the line of sight 
and the observed proper motions strongly support the opposite case. The first scenario points to a 
metal-poor lens with mass M » 0.6M Q at the distance D\ ss 4 kpc. In the second scenario the lens 
could be a main-sequence star with M = 0.8-0.9M Q about half-way to the Galactic bulge or in the 
foreground disk, depending on the extinction. 

Subject headings: gravitational lensing — Galaxy: center — Galaxy: bulge 
— stars: fundamental parameters (colors, masses) 



1. INTRODUCTION 

Gravitational m icrolensing of st ars within the Local 
Group of galaxies (jPaczvhski II 1996( 1 directly probes both 
luminous and dark matter concentrations along the line 
of sight. Over the past decade microlensing surveys 
have continued to enable observations with far-reaching 
implications, such as constraints on the fraction and 
content of Galactic dark matt er (e.g. lAlcock et al. I 
119961: lAlcock et af] 119981 120011 ). di scovery and char- 
acterization of exo - planet systems dBond et al.l 120041 
Udalski et ail l2005t iBeaulieu et"aH 120061 : iGould et all 
20061) . and measurements of the fundame ntal properties 
of stars and their evolutionary end point s (|Bennett et al.l 
120021 1 Abe etaD200alGould 'et al.ll2004f ). Unfortunately, 
while the light curve of a microlensing event provides the 
key discovery signature, it is insufficient to solve uniquely 
for the mass, the distance and the relative transverse 
velocity of the lens. As a result, out of a few thou- 
sand events discovered to date, only a handful a l lowed 
the mass of the lens to be measur ed (|An et al. I l2002t 
IGould et al.ll200ll.Iiang eTaD 12003 ) . 

In the case of microlensing by a luminous body (a star) 
the basic degeneracy of the model can be broken by di- 
rectly observing both the lens and the source. The diffi- 
culty with this approach, however, is inherent in the ge- 
ometry of microlensing that implies milli-arcsecond sep- 
arations between the lens and source components dur- 
ing the event. So far MACHO-LMC-5 was the only mi- 
crolensi ng event for which the lensing body has been re- 
solved (|Alcock et al.l l200l| ) . The lens that gravitation- 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele- 
scope Science Institute. STScI is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. 

2 Jodrell Bank Observatory, University of Manchester, Maccles- 
field, Cheshire SK11 9DL, UK, e-mail: simkoz@mlens.net 

3 Los Alamos National Laboratory, MS-D466, Los Alamos, NM 
87545 



ally magnified the source in the Large Magellanic Cloud 
turned out to be a nearby M dwarf in t he Galactic disk 
([Drake et all 12004 IGould et all 120041 ). iBennett et aT~l 
(2006) demonstrated the presence of a bright lens com- 
ponent in the planetary microlensing event OGLE-2003- 
BLG-235/MOA-2003-BLG-53 and estimated the mass of 
the host star using the centroid shift of the combined 
light. 

Here we report a direct detection and mass mea- 
surement of the gravitational lens responsible for the 
MACHO-95-BUL-37 event— the first fully resolved mi- 
crolensing system involving a Galactic bulge source, and 
the second such system in general. 

2. MICROLENSING EVENT MACHO-95-BLG-37 

The event was discovered by the MACHO collabo- 
ration as a single and apparently achromatic bright- 
ening of object 109.20635.2193 in their photomet- 
ric monitorin g dat abase of the Galactic bulge (see 
iThomas et"aD 120051 ). The object is quite faint (V ~ 
20 mag) and located in one of the densest fields cov- 
ered by the survey: equatorial (J2000) and Galac- 
tic coordinates (a, (5) = (18h04m34.44s, -28°25'33.7"), 
(/, b)=(2°54, -3°33). The location of the MACHO-95- 
BLG-37 event was one of the targets in our proper mo- 
tion mini-survey of the Galactic bulge ({Kozlowsk i et al.l 
120061 ) . Each of the 35 fields in the mini-survey was cen- 
tered on a microlensed source and was covered by two 
Hubble Space Telescope (HST) pointings taken several 
years apart. Using several relatively isolated stars we 
could co-register the HST and ground-based MACHO 
images to within 0.1" and unambiguously identify mi- 
crolensed sources, even in the presence of additional stars 
that were only resolved in the HST images. 

In the case of MACHO-95-BLG-37 we found that not 
only is the microlensed source accompanied by another 
very close star with comparable brightness, but also that 
the relative proper motion of the two components places 
them within 2.6 ± 3.5 mas of each other on 21 September 
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Fig. 1.— HST images of the MACHO-95-BLG-37 event location. 
The first epoch WFPC2/PC image (left) and the second epoch 
ACS/HRC image (right) were taken, respectively, 3.71 and 8.95 yr 
after the maximum light. The relative trajectory of stars A and 
B places them within 2.6 mas of each other at the time of the 
microlensing event. The circle shows a seeing disk characteristic of 
the ground-based MACHO survey (FWHM~ 2".l). The crosshair 
indicates the unbiased centroid of the lensed light from difference 
imaging. 



1995 (HJD 2449982.3) when the microlensing event took 
place. The prior probability that the blend is a random 
coincidence is very small, so we have a clear indication 
that the companion source is actually the gravitational 
lens that caused the event of 1995. Before we begin a 
detailed investigation of this finding (§§[3]-[5j) we first de- 
scribe the available data and basic data reductions. 

2.1. HST astrometry and photometry 

A detailed description of the rele vant HST data 4 was 
published in iKozlowski et alj (|2006l ) and only the essen- 
tial facts are repeated here. The first and second epoch 
images were collected, respectively, 3.71 and 8.95 yr af- 
ter the event. The first pointing employed the Plane- 
tary Chip (PC) of the Wide Field Planetary Camera 
2 (WFPC2) instrument, and provided (nearly) simul- 
taneous color information in both V and I photomet- 
ric bands using F555W and F814W filters. During the 
second pointing we used the High Resolution Channel 
(HRC) of the Advanced Camera for Surveys (ACS) and 
obtained high signal-to- noise ratio (S/N) imaging in the 
F814W filter only. In each case we co-added all suitable 
F555W and F814W images for a given epoch. The field 
of view covered by the ACS/HRC and WFPC2/PC de- 
tectors is similar (29" x 26" and 35" x 35", respectively), 
but not the pixel size (25 versus 45.5 mas). At the HST 
resolution the MACHO database object associated with 
the microlensing source was immediately revealed to be 
a composite of four unresolved stars, which we label A 
through D (Fig. [TJ. 

The magnitudes and positions of stars A-D were ex- 
tracted from the fits of stellar profiles. The local point 
spread function (PS F) models were ge nerated using the 
TINYTIM software (jKrist 1119931 119951 ) and interpolated 
with bi-cubic splines. For all model fitting we used the 
MINUIT package 5 . Stars A and B have overlapping pro- 
files and required a special model with two PSF com- 
ponents fitted simultaneously. A small section of the 
ACS/HRC image was fitted first, providing an unbiased 
value of the second epoch separation between the two 
components and a good handle on the flux ratio in the 

4 http: //archive . stsci ■ edu/h stT] 

5 http: / /wwwasd. web . cern. ch/wwwasd/cernlib/ 
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Fig. 2. — Light curve of the MACHO-95-BLG-37 event from our 
Difference Image Analysis of data in the public MACHO database. 
A significant color change near the peak magnification indicates 
that the microlensed source is redder than the unresolved compos- 
ite of stars A, B, C and D from Fig.fi] 



^F8i4W _ band. The flux ratio was then fixed at the sec- 
ond epoch value for the purpose of fitting the /f814W _ 
band WFPC2/PC image and obtaining the first epoch 
astrometry. Finally, the (V — I) colors of stars A and 
B were established by fitting the VF555W _ band image us- 
ing a model with variable flux ratio and the blend sep- 
aration fixed at the value taken from the /F8i4W _ band 
fit for the same epoch. The resulting astrometric and 
photometric measurements are given in Tables [1] and [H 
Note that in the F-band the only available high resolu- 
tion imaging comes from the relatively shallow first epoch 
WFPC2/PC observation, so the A/B flux ratio is poorly 
constrained and the errors in V and (V — I) are rela- 
tively large for these two stars. In the /-band, however, 
we have an accurate measurement of the flux ratio from 
ACS/HRC that allowed us to eliminate a degenerate free 
parameter from the double PSF fit of the WFPC2/PC 
data. This explains why the astrometric accuracy for the 
first epoch is actually better than for the second epoch, 
despite a larger pixel size and a much smaller separation 
between stars A and B in the WFPC2/PC images com- 
pared to the ACS/HRC data. Using simulated images 
we found that our first epoch astrometry is biased by 
about 1.5% toward lower separations. We could not find 
a better procedure that would eliminate this effect, so a 
post-factum correction was included in the WFPC2/PC 
data reported in Table [1] 

2.2. Microlensing light curve revisited 

The MACHO-95-BLG-37 event was recorded in a faint 
star subject to intense crowding, and therefore the stan- 
dard light curve in the MACHO photometric database 
has a very low S/N. In order to reduce the uncertainties 
of the microlensing parameters derived from light curve 
modeling we performed Difference Image Analysis (D IA ; 
lAlard fc Luptoril [j~99ll lAlardll^pa IWozniak 112000) on 
the original ground-based images, i.e. on the simulta- 
neous two-color imaging data collected by the MACHO 
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Instrument Epoch t — to Aa AS 

(yr) (yr) (mas) (mas) 

WFPC2/PC 1999.43 3.71 -29.9 ±1.3 -24.9 ±1.3 
ACS/HRC 2004.67 8.95 -74.5 ±1.6 -63.1 ± 1.6 

Note. — Aa, A<5 are positions of star B relative to star A 
(Fig. [T} in a Cartesian reference frame aligned with the local equa- 
torial coordinates. The moment of maximum light to corresponds 
to HJD = 2449982.3 (Epoch 1995.72). 



Table 2. HST photometry 



Star 


VF555W 
(mag) 


-7F814W 
(mag) 


(V-I) 
(mag) 


fv 


fi 


A 


20.24 ±0.08 


18.45 ±0.02 


1.79 


0.28 


0.35 


13 


20.30 ±0.08 


19.07 ± 0.04 


1.23 


0.26 


0.19 


C 


20.29 ±0.03 


18.66 ±0.03 


1.63 


0.27 


0.28 


D 


20.63 ±0.04 


19.18 ±0.04 


1.45 


0.19 


0.18 



Note. — fv, fi are fractional contributions to the total flux 



survey 6 . The PSF matching and photometric solutions 
were confined to a 2'.5 x 2'.5 region around the source 
(approximately 256 x 256 pixels). After discarding ob- 
servations outside the relevant time interval and rejecting 
a small fraction of frames with bad seeing we considered 
a total of 132 images in each of the MACHO photometric 
bands 6m and tm- High S/N reference images were con- 
structed by co-adding 9 good quality images with a well- 
behaved PSF. From a series of difference frames in which 
the source was significantly magnified we derived an un- 
biased centroid of the lensed light that clearly points 
to the pair of stars A and B when transformed to the 
ACS/HRC coordinates (Fig. [1]). One of these two stars 
must then be the microlensed source. 

The reference flux in each band was derived from a 
comparison between our differential fluxes and conven- 
tional PSF photometry obt ained with the Dophot soft- 
ware (jSchechter et al.l ll993) running in a fixed-position 
mode with the input object lists based on our deep ref- 
erence images. Stars C and D could not be properly 
deblended, even using fixed HST positions transformed 
to the template coordinates. The template position of 
the A-D composite was set to the mean ACS /HRC posi- 
tion of stars A and B. We selected 31 calibration images 
per photometric band with the best overall seeing, back- 
ground and transparency. In seven of these images the 
source was visibly magnified. The statistical uncertainty 
of the reference flux is 8% in 6m and 9% in tm- The 
background level estimated by the Dophot algorithm in 
a crowded field is somewhat sensitive to the assumed 
shape of the PSF (especially in the wings). In our case 
of a very faint object near the detection limit set by the 
source confusion we find that the systematic uncertainty 
in the reference flux can easily reach 10%. This generic 
problem is partially alleviated by the fact that the sys- 
tematics are similar in both filters and source blending 
must always be considered in the analysis of individual 
light curves in crowded fields. 

6 http: //wwwmacho .monaster . ca/Data/MachoData. html 



Parameter Value Error 



t (days) 982.3 0.3 

t E (days) 25.2 4.2 

u 0.37 0.10 

f aV 0.33 0.12 

f s ',R 0.38 0.14 

m v (mag) 19.314 0.005 

m R (mag) 18.545 0.003 

xl 1-490 

v 255 



Note. — Maximum magnification is at to days after HJD = 
2449000. 



The final light curves (Fig. ^ were shifted to the in- 
strumental 6m, ?"m scale of the MACHO database using a 
median offset for a few tens of bright stars near the loca- 
tion of the MACHO-95-BLG-37 event and transformed 
to approximately stan dard Vm , Rm magnitudes following 
iPopowski et al.l (|2005f) . We also determined transforma- 
tions between Vm, Rm and the standard V, I magnitudes: 

V M = V + (0.05 ±0.11)(V-J) + (0.01 ±0.20), m 
R M = I + (0.62 ± 0.10) (V — I) — (0.17 ± 0.18). [i) 

Hereafter, the subscript is omitted and MACHO filters 
are implied for V, R photometry. The overall quality of 
the U-band light curve is lower compared to that in the 
i?-band due to occasional pixel level defects in the 6m 
frames that were clearly visible in the difference images. 

3. MICROLENSING LIGHT CURVE MODELS 

The first step is to obtain the basic microlensing pa- 
rameters such as the time-scale tg, the dimensionless im- 
pact parameter uq, the moment of the peak brightness 
to, and the baseline magnitudes my.R- hi order to pre- 
serve consistent color information, both V- and i?-band 
light curves were fitted simultaneously with a simple mi- 
crolensing model that allows for flux blending (source 
fractions / s < 1). The data point at t = 1016.9 days 
is a moderate outlier in the U-band light curve (Fig. \2§ 
and is rejected in all analyses. The change in \ 2 due to 
this cosmetic change is not significant and none of our 
conclusions are affected. The resulting best fit model is 
given in Tableland provides a marginally acceptable fit 
(reduced \1 — 1-49 for v — 255 degrees of freedom). 

3.1. Colors 

Using different parameterizations of the model equiva- 
lent to the one in Tabled we obtained the source/blend 
colors and the color difference with the error bounds that 
fully account for covariance: (V — R) s — 0.92 ± 0.04, 
(V - i?) b = 0.68 ± 0.04 and A(V - R) s , h = (V - R) s - 
(V — i?,)b = 0.24 ± 0.06. This corresponds to a positive 
color shift during the event A(V — R) cvcn t — ±0.06 mag 
and indicates that the source is redder than the blend. 
However, it must be emphasized that the measurement of 
the reference flux for our light curves poses a significant 
challenge given the limitations of the available archival 
data (§EH). Both (V - R) h and A(V- R) s , h are subject 
to the systematics of the reference flux in two bands. The 
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value of (V — R) s , on the other hand, is more reliable, 
because it is constrained by the magnified portion of the 
light curve, even if the reference fluxes are not known. 
This is best seen from the model of the simultaneous two- 
color DIA light curve written as AF(t) — F s x A(t) + F a 
in each band, where A(t) is the magnification factor and 
Fq < if the source is effectively magnified in the refer- 
ence image. Although in most cases the source flux F s 
is poorly constrained in both colors, the error bounds on 
the ratio F B v/F Bt R are relatively tight due to covariance 
and, most importantly, independent of the flux offsets. 
Therefore, the derived value of (V — R) s only depends on 
the global calibration of flux units for the reference im- 
ages, which can be done much more reliably using bright 
isolated stars. In conjunction with the HST photometry, 
the source color information will be crucial to deciding 
the identity of the microlensed source, and therefore the 
lens (S OI) . 

3.2. Parallax constraints 

The ground-based microlensing light curve provides 
useful constraints on the acceleration term in the ob- 
served trajectory of the lens relative to the source. In 
the case of a short, low- magnification microlensing event 
such as MACHO-95-BL G-37 we can only obtain one- 
dimensional information (iGould et al.lll994l ). Following 
the geocentric formalism of IGould I (I2004TT we introduce 
into the model the dimensionless microlensing parallax 
vector 7Te, where 7Te,|| is the component of 7Te opposite 
the direction of the projected position of the Sun at the 
peak of the event. We find that 7r E .|| = 0.07tQ4g, while 
tte,j_ remains unconstrained, i.e. there is no detectable 
parallax. There are two observations during the event 
(at t = 967.9 and t = 996.9 days) with atypically low 
y-band fluxes and relatively large error bars compared 
to the adjacent measurements. Without these two data 



points we get 7Te 



0.00 



+0.67 
-0.45 



and the apparent weak 



asymmetry of the best fit model goes away. In ij l4.ll this 
constraint is improved using HST photometry and in 
§ 14.21 combined with the HST astrometry to place the 
limits on the relative source-lens parallax 7r re i. 

4. RESOLUTION OF THE MICROLENSING SYSTEM INTO 
LENS AND SOURCE 

The fundamental difficulty with resolving a lens de- 
tected through time- variable magnification is that its ap- 
parent separation from the source is below the HST res- 
olution for months or even years after the event. In the 
case of MACHO-95-BLG-37 (and similarly for MACHO- 
LMC-5) this problem is greatly reduced due to the rather 
large relative motion of stars A and B (§[13}. High pre- 
cision HST astrometry at two epochs well after the peak 
magnification allowed us to calculate a very accurate rel- 
ative trajectory of star B with respect to star A. Simply 
connecting the two measurements in Table Q] we get: 



r(t) - r(h) 



(*2 - t) 



(t-ti) 



(2) 



(t 2 -ii) v y (* 2 -ti)' 

where r = (Aa, AS) is the relative position with mea- 
surements available at t\ — 3.71 and £2 = 8.95 yr. The 
separation at the peak of the event (t = 0) is then: 

(Aa , AS a ) = (1.6 ± 2.5, 2.1 ± 2.5) mas. 
These values are fully consistent with a model in which 
the two stars are the source and lens, and which pre- 



dicts a very low value of the two-dimensional separation 
ro = r(0). There are two alternative possibilities: cither 
one of the members of the pair is a random interloper, 
or it is a companion to either the lens or the source. The 
first possibility is ruled out by the following argument: In 
the sky region under consideration the density of stars is 
0.085 and 0.176 per square arcsecond for stars brighter 
than I = 18.45 and 19.07 mag, respectively. The cor- 
responding Poisson probabilities of a random alignment 
within 2.6 mas at the time of the event are 1.8 x 10~ 6 
and 3.7 x 10 -6 , respectively, i.e. very low. The other 
case, of one of the two detected stars being a compan- 
ion to either the lens or the source, can also be ruled 
out. It is clear is that one of the two stars must be the 
source. Furthermore, the rapid relative proper motion 
excludes the possibility that the second star is a com- 
panion of the source (the implied binary motion will be 
too high, about 400 km/s at a distance of 8kpc). Thus 
we only need to consider the possibility that one of the 
stars is a companion of an unseen dark lens, with a sep- 
aration of about 2.6 mas between them (recall that the 
lens is almost perfectly aligned with the lensed source 
at the peak). This is about b w 3.5 Einstein radii (for 
#e = 0.75 ± 0.13 mas, see §4.2). We can approximately 
model the perturbation of the luminous star on the dark 
lens as a Chang- Refsdal lens (Chang & Refsdal 1984). 
The shear induced by the luminous star at the position 
of the dark lens would be 7 = q/b 2 » 0.08, where q is 
the mass ratio of the luminous companion to the dark 
lens. The short t~E does not favor a massive dark lens 
such as a black hole and neutron star, and so the mass 
ratio q is likely larger than one. The caustics will have a 
size roughly « 2q/b 2 ~ 0.16 (e.g. Mao 1992). 

The caustics size is comparable to the measured impact 
parameter (uq ~ 0.37), which would introduce a strong 
asymmetry in the light curve for most trajectories (not 
seen in the observed low S/N light curves). Hence we 
regard the 'dark' lens scenario as not very likely. The 
bright lens hypothesis is thus favored and we conclude 
that the source and lens system involved in the MACHO- 
95-BLG-37 microlensing event consists of stars A and B 
from Figure [1] (in an order still to be determined) . Our 
subsequent arguments are based on that assumption. 

4.1. Identifying the lens 

To find out which member of the candidate pair of stars 
is the lens, we can make use of the observed color change 
during the event and the fact that gravitational lensing 
is achromatic. In Figure [3] we plot the color-magnitude 
diagram (CMD) of the stellar field around MACHO-95- 
BLG-37 and stars A-D from Figure [TJ The light curve 
in Figure [2] reflects the integrated flux of the four stars 
(unresolved in ground-based images). From the observed 
color increase A(V — R) e vent — +0.06 mag near the peak 
magnification we infer that the microlensed source is red- 
der than the composite. Although the (V—R) colors were 
not measured individually for stars A-D, it is very un- 
likely that the ordering of the (V — I) and (V — R) colors 
is different. Stars A and B are, respectively, the reddest 
and the bluest components of the blend, so the color shift 
points to star A as the source. This is entirely in agree- 
ment with the source color (V — R) s = 0.92±0.04 mag 
and the color difference (V — R) s ,h — 0.24 ± 0.06 mag be- 
tween the source and the rest of the blend found in §13.11 
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Equation fl} implies HST (V — I) s ~ 1.7 mag, also con- 
sistent with the HST photometry of star A. Thus, based 
on the light curve evidence, star B must be the lens, be- 
cause star A is the microlensed source. However, in §[5] 
we show that the physical interpretation of the CMD and 
kinematic data strongly argues against this result. 

In principle, the HST photometry provides an addi- 
tional test of these possibilities because we can trans- 
form the measurements to the MACHO system and 
obtain a constraint on the source magnitudes V = 
20.34 ± 0.09, R = 19.38 ± 0.05 assuming star A and 
V = 20.37 ± 0.11, R = 19.66 ± 0.08 assuming star B 
(including the variance and covariance in transformation 
coefficients). Unfortunately the difference A% 2 ~ 0.1 
between the models with stars A and B as the source 
is completely insignificant. An additional problem is 
that the contributions of stars A-D to the blend are not 
known very well. In order to match the total MACHO 
baseline magnitudes (§ 12. 2| we would have to make the 
transformed HST fluxes of all four stars fainter by 10- 
15%, depending on the photometric band, and then still 
assume that only about half of the flux in stars C and 
D is effectively added to the total flux of stars A and 
B. This is not surprising knowing that stars C and D 
are near the edge of the FWHM disk of stars A and 
B (Fig. [T]) — yet they are too faint to be deblended — and 
that in ground-based microlensing images the "sky" level 
is set by a featureless continuum of merging stars. Over- 
estimating the background by a mere few counts makes 
stars near the detection limit appear 0.1-0.2 mag fainter. 
The weights are probably slightly different in each pho- 
tometric band due to details such as the orientation of 
stars C and D with respect to the PSF that is never per- 
fectly round. Nevertheless, it is still useful to perform 
a microlensing light curve fit with a single additional 
"measurement" of the source magnitude, i.e. effectively 
constrain / s . Compared to the results in § 13.21 the er- 
ror bounds on the dimensionless parallax are improved, 
yielding 7r E || = 0.0±0.4 for any reasonable set of flux 
weights. 

4.2. Proper motion, Einstein ring radius and relative 

parallax 

The relative proper motion /x re ] of the lens with respect 
to the source provides further clues about the physics 
of the event. For the case at hand, a robust value of 
^ ro i can be derived from the HST astrometry alone. If 
we ignore the parallax and approximate uq — 0, we 
find ^ ro i = 10.85 ± 0.16 mas yr -1 and the position an- 
gle (j) = 229.83 ± 0.84 deg (North through East). Al- 
lowing for a finite impact parameter makes no differ- 
ence to /i re i, changes by a mere 0.4 deg, and pre- 
dicts the peak brightness a couple of months after the 
actual event, consistent within la uncertainties. Includ- 
ing the parallax also has a negligible influence on the 
trajectory. Thus the Einstein radius can be estimated as 
9e = Mrei*E = 0.75 ± 0.13 mas. 

The direction of the dimensionless parallax vector tte 
is the same as the direction of the lens-source relative 
proper motion. The component 7r E || that points away 
from the projected position of the Sun is almost perfectly 
due East, since the event peaked on 21 September. Using 
results from § 13.21 and § !4.1l we can immediately estimate 
= |tt e ii (sin 0) 1 1 < 0.53, and set an upper limit on 
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Fig. 3. — Color-magnitude diagram (CMD) of stars around 
the MACHO-95-BLG-37 event location. The HST photometry of 
stars A-D from Fig. [T] (open circles) is plotted against the average 
bulge and disk pop ulation along the line of sight (OGLE-II data; 
lUdalski et"aTll200l ). 

Table 4. HST kinematics of stars A and B 



star 


Ml 


Mb 




( mas yr — - 1 ) 


( mas yr — 1 ) 


A 


6.9 ±0.3 


-0.2 ±0.3 


B 


-3.8 ±0.4 


3.7 ±0.4 



Note. — Proper motio ns are expressed in an average star ref- 
erence frame as defined in Kozlowski ct al. (2006). 



the relative lens-source parallax, 7r rc i = 7Te#e ^5 0-3 mas. 

5. MICROLENSING SCENARIOS AND THE LENS MASS 

The value of the source color derived in § 14.11 favors 
a scenario in which star A is the source and star B is 
the lens. However, as we show in this section, such ar- 
rangement is very unlikely in the context of the CMD 
(Fig. [3]) and proper motions measured relative to the 
Galactic bulge (Table [4]) . Although the extinction-to- 
reddening ratio in the direction of the event is abnormally 
low, t he reddening anomaly cannot explain the conflict. 
After ISumi I (|2004f ) , we adopt the reddening coefficient 
i?vi = -A-y/Ey-i — 1-98 and the total bulge extinction 
Ay = 1.54 mag. In the following discussion we consider 
both source star cases in some detail and then use the 
measurement of #e from § 14.21 to constrain the mass of 
the lens. 

5.1. Blue lens scenario 

First we attempt to reconcile all available data with 
the evidence in § !4.1l that the source is red. Given its red 
color, star A is too faint to be a giant and too bright 
to be on the main sequence in the Galactic bulge. If 
it were a giant several magnitudes behind the bulge, it 
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could in principle belong to the Sgr dwarf galaxy, but 
its observed prope r motion is not consistent with Sgr 
(jlbata et al. 1119 97). More likely, star A is a dwarf in the 
foreground disk at a distance of ~ 2.5 kpc and behind 
most of the extinction. Then, if star B is indeed the lens, 
it must be in front of star A, and the only simple solution 
is that the lens is a nearby white dwarf at ~ 100 pc or 
so. Unfortunately, this exciting possibility is ruled out 
by the parallax constraint 7r ro i < 0.3 mas (§ 14.20 . as it 
predicts 7r rc i ~ 10 mas. 

The location of star A in the CMD is still marginally 
consistent with a faint subgiant on the far side of the 
bulge subject to ~ 0.2 mag of extra reddening compared 
to the general population. But there is little support for 
that, since the CMD shows a compact red clump and 
indicates a very unifor m extinction across this field (c.f. 
the extinction map of lSumill200l . The observed kine- 
matics would also be very unusual for this scenario with 
star A showing a 7 masyr -1 disk-like prograde motion 
in the plane and star B moving at a ~ 135° inclination. 
One could still argue that star B is a low metallicity halo 
subdwarf to explain its motion and dramatically increase 
the Ay prediction, but there is simply too much fine tun- 
ing to consider this a reliable solution. 

5.2. Red lens scenario 

The properties of both stars are much easier to explain 
if we dismiss for a moment the source color evidence from 
§ 14.11 and assume that star B is the source and star A is 
the lens. In this case star B is most likely in the bulge, 
where its absolute magnitude and color would be approx- 
imately Mi — 3.6, (V — I)o = 0.4 assuming a red clump 
at Jrc = 15.3, (V — J)rc = 1-8 in Fi gure [Hand ado pting 
Mi jRC = -0.2, (V - I)o,rc = 1.0 l)Udalski|[200l . So 
the source fits the properties of a metal-poor star near 
the turnoff point in the bulge, and the observed proper 
motion is fully consistent with this picture. Then star A 
must be the lens and can be placed on the main sequence 
at a distance of ~ 4 kpc, where it would follow the Galac- 
tic rotation near the plane and move a few masyr -1 . 
Again, the observed kinematics support this scenario. 

5.3. The lens mass estimates 

Any acceptable scenario for the lens must satisfy the 
constraint on microlcnsing geometry set by the measure- 
ment of the Einstein radius 6% = 0.75 ± 0.13 mas (§ 
For the lens of mass M we have 



M = 



Mr. 



3.14 mas, 



K7T rc l 



(3) 



where 7r 



rcl 



d; 



D s 1 is the relative parallax for the 



lens and source distances D\, and D s kpc. A given value 
of the source distance D s sets a relationship between the 
lens mass M and the range of lens distances allowed by 
the error bounds of #e- Making a reasonable assumption 
about the luminosity class of the lens we can parameter- 
ize the photometric solutions in the same way, i.e. using 
the mass of the lens. 

For each value of the lens mass M we use the appro- 
priate mass-luminosity-color relation to obtain the abso- 
lute magnitude My and color (V — I)q. Then using the 
HS T photometry (§ 12. 1[) we estimate the reddening -EV-i 
through Ey-i — (V— I)hst — (V— I)o, and extinction Ay 




Lens A 

- — Lens B 
— ■ Lens B [Fe/H] = -1 



0.4 



0.6 



0.8 
M/IV1 

Sun 



1.2 



Fig. 4. — Photometric solutions for the lens distance and required 
extinction (top panel) as a function of the lens mass are shown. 
The gray area covers microlensing geometries for a source at 8 kpc 
allowed by the measurement of the Einstein radius #e = 0.8 mas 
with 20% uncertainty. For each main-sequence model we show a 
pair of lines that reflects the ±0.1 mag uncertainty in stellar colors. 
A metal-poor model is also shown for star B (dot-dashed line) . The 
top panel shows the required extinction to match the observed 
magnitude and color of the lens f £| 5 - 3 I t ; the thick horizontal line 
marks the to tal Galactic e xtinction (Ay = 1.54 mag) for the line 
of sight from lSuml 



using the reddening coefficient Ry i = Ay/Ey i = 1.98 
taken from the extinction map of lSumil d200l . Com- 
bined with the assumption of the source located at 8 
kpc, each set of the above parameters allows a calcu- 
lation of the lens distance D\ and extinction Ay, which 
should not exceed total extinction of the bulge Ay = 1.54 
mag (see Fig. . We adopted th e mass-luminosity rela - 
tion for the main sequence from ISchmidt-Kaler I (fl982T ) 
and the empirical color-magnitude relation defined by 
our polynomial fit to the Hipparcos C MD data in abso- 
lute m agnitudes (Hipp arcos catalogue; iPerrvman et al.l 
(1997), iBessell I (fl990ft . compiled by I.N. Reid 7 ). Our 
CMD locus for the main seq uence is very close to the 
linear relation of iReid I (|1991| i for O.5M < M < 1.0M Q 
and is brighter by up to 0.3-0.5 mag outside this range. 
For comparison we also used a model grid of low metal- 
licity hydrogen-bur ning stars with [Fe/H] = —1.0 from 
iBaraffe et all (|1997l ). 

The constraints on the lens resulting from the two ma- 



http: //www-int . stsci . edu/^inr/cmd.html 



First direct /u-lens toward the bulge 
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jor scenarios are plotted in Figure [4j The red solid lines 
show the photometric constraints for star A being the 
lens and the blue dashed lines are for star B being the 
lens. Each pair of lines corresponds to a range of so- 
lutions allowed by the ±0.1 mag uncertainty in stellar 
colors. The gray area is the geometric constraint based 
on the measurement of #e (eq. [5]) assuming a Galac- 
tic bulge source. It is clear that a blue lens (star B) on 
the main sequence generally underpredicts the amount 
of extinction for a given distance. If the source is in the 
bulge, the microlensing constraint selects Dj=4-5 kpc, 
where the blue color of the lens does not allow any signif- 
icant extinction. Moving the source to a distance much 
larger than D s — 8 kpc shifts the range of D\ upward 
by a few kpc, but the problem with Ay remains. The 
blue dot-dashed line illustrates the effect of lowering the 
metallicity of star B to [Fe/H] = —1.0. This does solve 
the issue of extinction, but requires that the lens is a 
metal-poor subdwarf in the Galactic halo, or perhaps in 
the thick disk. Such a possibility is unlikely since the 
halo and thick disk contribute only a small fraction of 
stars within the Galactic disk. The mass and the dis- 
tance of the lens are then M ~ 0.6M Q and D\ ~ 4 kpc. 
The solution with a red lens (star A) is also not without 
a wrinkle, because in order to avoid overshooting the to- 
tal extinction for the bulge we need to make #e about 
20% larger and the lens ~ 0.1 mag bluer compared to 
the best estimates. Nevertheless, we can still find a con- 
sistent answer within la uncertainties. In this scenario 
the lens is a main-sequence star with M = 0.8-0.9M Q 
(spectral type G5-K0) at a distance of D\ ss 4 kpc. 

6. DISCUSSION 

There is little doubt that we are directly observing 
the lens in the MACHO-95-BLG-37 event as it separates 
from a nearly perfect alignment with the microlensed 
source. However, the final identification of the gravi- 
tational lens is somewhat problematic. While the light 
curve models combined with the photometric data for in- 
dividual objects favor a scenario with a blue lens and a 
red source, the opposite assignment is much more plau- 
sible in the context of the color-magnitude diagram for 



the line of sight and the observed proper motions. In any 
case, the lens is a relatively bright star with the mass of 
~ 0.6M© or ~ 0.9M Q . It is conceivable that additional 
factors such as binarity of stars affect the interpretation 
of the MACHO-95-BLG-37 event, but a conclusive reso- 
lution of the present conflict will require new data. High 
quality spectroscopy would unambiguously pinpoint both 
the 3-D kinematics and the distance scale. 

For the first time in a Galactic bulge event the lens 
and source have been directly resolved. This is an im- 
portant addition to the sample of one consisting of the 
MACHO-LMC-5 event. There are thousands of known 
Galactic bulge microlensing events and several dozen of 
those have archival HST pointings suitable for follow-up 
proper motion work. I n our proper motion mini-survey 
(Kozlows ki et al.l l2006| ) we included 35 of those fields, 
and have already identified sever al more promising can - 
did ate lenses. A s pointed out by lHan &z Chang 1 (2003) 
an dWood (2006), in a few per cent of microlensing events 
toward the Galactic center a lens with characteristic mo- 
tion /z rc i < lOmasyr -1 may be detectable a decade af- 
ter the microlensing episode. We are only beginning to 
probe directly the mass spectrum of the Galactic mi- 
crolenses; however, we can expect that in the short term 
the progress will accelerate considerably due to availabil- 
ity of the archival and future HST pointings. 
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